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Abstract. We investigate the scaling properties of the excitation energies and transition ampli- 
tudes of the one- dimensional spin-^ antiferromagnetic Heisenberg model exposed to an external 
perturbation. Two types of perturbations are discussed in detail: a staggered field and a dimerized 
field. 

PACS. 75.10 -b General theory and models of magnetic ordering 



1 Introduction 

The one-dimensional spin-^ antiferromagnetic Heisenberg 
model with nearest neighbour couplings and periodic bound- 
ary conditions (Sn+i = Si): 



coupled spin chains. Treating the interchain coupling in 
the mean field approach, reduces the system to a one- 
dimensional Heisenberg model in a staggered field pfjj . 

It is the purpose of this paper to study excitation en- 
ergies 



N 



H 



(1.1) 



has been studied intensively with analytic and numerical 
methods. Eigenvalues as well as transition matrix elements 
for the spin operator have been calculated with the Bethe 
ansatz and quantum group symmetries These calcu- 
lations allow to exploit part of the dynamical properties 
of the model, the two-spinon contributions, at T = 
The dynamics of the model in the presence of an external 
magnetic field, which is periodic in space 



H(/i) =H + /iS 3 (g) 



with 



Sa (5) 



N 

i=i 



1 ' 



(1.2) 



(1.3) 



has been studied so far mainly for p = jM0,pLp],|^,||,p|]. 
Since the total spin 83(0) commutes with H, the eigenvec- 
tors of H(/i) and H are the same, the eigenvalues change 
in a trivial manner . Fo r this reason the magnetic proper- 
ties of the model (1.1) in a constant field - as there are 
magnetization curves, susceptibilities, static and dynamic 
structure factors - were accessible by means of the B ethe 
ansatz as well. The case of a staggered field, i.e. Eq. (1.2) 
for p = 7r, in a one-dimensional Heisenberg model has 
attracted considerable interest as an effective model for 



u mn (N, h) = E m (N, h) - E n (N, h), 



and transition amplitudes 



T mn (N,h) ee (& m {h)\S 3 {ir)\& n {h)), 



(1.4) 



(1.5) 



for the operator 83(77) in the one-dimensional antiferro- 
magnetic Heisenberg model in a staggered field of strength 
h. 

In section || we derive a system of differential equa- 
tions, which describes the evolution in h for u> mn and T mn . 
Section || is devoted to a finite-size scaling analysis of the 
ground state excitations uj m o an d T m Q in the combined 
limit h — * and N — > 00. It is shown in section ^, that 
the critical exponents a and e in the scaling ansatz and the 
scaling variable as well can be computed by means of the 
evolution equations. For small values of the scaling vari- 
able the scaling function is also determined. In section |^ 
we analyse in the same way the dimer operator. 



2 Evolution equation for excitation energies 
and transition amplitudes 



Star ting from the eigenvalue equation of the Hamiltonian 



H(h)\¥ n (h)) =E n (N,h)\V n (h)), 



(2.1) 
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it is straightforward to derive the evolution equations for 
the energy eigenvalues and eigenstates: 



dE n (N, h) 
dh 



T nn (N, h), 



>•■<"» - - £ 



T mn {N, h) 
i(N,h) 



(2.2) 
(2.3) 



where we have used the definitions ( |P| ) and (1.5). 

The evolution is governed by the tran sition matrix el- 
ements of the operator 83(77) [Eq. (L3)]. The operator 
83(77) is real hermitian and therefore the matrix (1.5) is 
symmetric T mn = T nm . 

Remarkably enough, there exists a closed system of 
differential equations, which describes the /i-dependence of 
the energy eigenvalues E n — E n (N, h) and the transition 
matrix elements T mn — T mn (N, h): 



d 2 E n 
dh 2 

dh 



-E 



IT, 



TmlTi n 



(2.4) 



Win 



Tmn dtO nrr 

dh 



(2.5) 



The Hamiltonian (1.2) behaves under translations T of 
one lattice spacing as: 



TH(/i)T f = H(-/i). 
Since T is unitary, the eigenvalues 

E n {N,h)=E n (N,-h) 
are symmetric, whereas the eigenstates |#(h)), 
\¥{-h))=T\9(h)), 



(2.6) 
(2.7) 
(2.8) 



are no longer momentum eigenstates, but linear combina- 
tions of two momentum eigenstates with momenta p = 
and p = 7T, respectively. Of course, in the limit h — > 
translation invariance is recovered and p is again a good 
quantum number. We use the following notation for the 
energy eigenstates: 



T\W n (0)) = ±W n (0)), 



0,1,2,. 



(2.9) 



which means these states carry momentum p — or p = 
7T, respectively. The ordering of the energy eigenstates is 
chosen in such a way that the ground state |#o(0)) belongs 

to the ordered sequence of eigenstates with n = 0, 2, 4, 

All eigenstates in this series have the same momentum. 
The series of eigenstates with n = 1, 3, 5, . . . carries the 
opposite sign. It starts with the first excitation ^(O)) and 
is independently ordered. 

Since the operator 83(77) changes the momentum of 
the state by 71", we get immediately from momentum con- 
servation: 



T m „(A^,0) = for (-1) 



m+n 



1. 



(2.10) 



The transition matrix elements T m Q{N, 0), m = 1, 3, 5, . . . 
from the ground state |5fo(0)) to the excited states |!^ m (0)) 
appear in the definition of the dynamic structure factor 

S(N,Lu,p = 7r,h = 0) = 

i£*[w-w^(JV,0)]|T^(^0)| 2 , (2.11) 

m 

which is known to develop an infrared singularity in the 
thermodynamic limit p[: 



D(00, Ul, 71", 0) — ► A . 



(2.12) 



The gap between the ground state and the excited states 
vanishes in the thermodynamic limit: 



Umo(N, 0) 



iV- 



+ 00 Q>mO 



(2.13) 



The coefficients a m o can be computed in principle by the 
methods of conformal invariance and finite-si ze sca ling [ O , 
O,0|. The next order terms entering in Eq. (2.13) are log- 
arithmic corrections. For this reason it becomes difficult to 
fix these coefficients from our system sizes up to N — 24. 
The difficulty of this task is extensiv ely discussed in ref- 
erences In figure £y(a) a fit without 
logarithmic corrections Jaw — 8.53) and with logarith- 
mic corrections (aio = 9.74) is shown. The exact value is 



9.86. 



o 





Fig. 2.1. Finite-size dependence of the ground state en- 
ergy (a) and the corresponding transition amplitude (b) of 
Hamiltonian (1.2), for system sizes N = 8, 10, . . . ,24 (•) The 



solid lines represent fits to the asymptotic behaviour (2.12) 
and (|2.14[). 



Neglecting the logarithmic corrections in the large- N 
behaviour, we expect for the transition amplitudes: 



T m0 (N, 0) N ^ b m0 N K , m = 1,3, 



(2.14) 



The exponent k turns out to be k ~ 1/2, as can be seen 
from the numerical data for Tiq(N, 0) in figure [2.l| (b). This 
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is in agreement with conformal field theory, from which 
in the present case follows K = am/27r 2 (Se e |19|,P0| and 
r efere nces therein). We consider Eqs. ( 2.10| ), J2.13| ) and 
( 2.14 ) as initial con ditio ns for the system of differential 
equations (^3) and (2.5). 



3 The finite-size scaling analysis 

In this section we present numerical evidence for the va- 
lidity of a finite-size scaling ansatz: 



w m o{N, h) 
uJ m a{N,0) 



1 + e m0 (x), 



T m0 (N, h) 



1 + fmoi.x), m=l,3, . 



T m0 (N,0) 
in the combined limit 

N — > oo, h — ► 0, with fixed x 



Nh e 



(3.1) 



(3.2) 



(3.3) 



Note, standard finite-size scaling theory 21,22] usually 
employ y = x 1 ^ as the scaling variable. For our purpose 
it is more convenient to use x instead. 

In figures |3.1| a nd |3.2| we show our numerical results 
for the ratios ( |3.1| ) and j3.2| ) on finite systems with N = 
8, 10, ...,24. The best results are achieved for 

e = 0.66(1). (3.4) 

The small z-behaviour of the scaling functions is well 



1.2 



o 
Z' 



1.1 



e 10 (x) = 0.0196(5) x 2 %(5) 


/ » 






e 20 (x) = 0.00 32(5) x 3 00(5) 








e 30 (x) = 0.0017(5) x M{S) y 
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Fig. 3.1. The ratio (fyj) for rn 
c 



1,2,, 6 versus {Nhy with 
0.66. The solid lines show the fits to the small x-behaviour 



of e m0 (x) [cf. (3.5) 



parametrised by: 



fmo{x) = fmOX* f , 



(3.5) 
(3.6) 



o 

g 

o 
E 

Z 



0.9 



0.8 



0.7 



f |0 (x)=-0.0177(15)x 



2.94(5) 



f, (x)=-0.0075(15)x 



2.99(5) 



10 



15 



(N h E ) 3 



Fig. 3.2. The ratio (fT|) for m = 1,3 versus (Nh e ) 3 with 
e — 0.66. The solid lines show the fits to the small a;-behaviour 
of /mo (as) [cf. (0)]. 



with 



4> e = 4>f = 3.00(4). 



(3.7) 



The diagonal matrix elements of the stag gered operator 
S 3 (tt) vanish in the limit h -> [cf. Eq. Q2.10fl]. There- 
fore, a scaling ansatz of the type ( |3.2D does not make 
sense. However, according to ( |2.2| ) the diagonal matrix el- 
ements of S3(7r) can be identified with the first derivative 
of the corresponding energy ei genv alue. Assuming, that 
these scale in the same manner (3.1) as the excitation en- 
ergies, we expect: 

-2 



Too (N, h) -> Nh° f 00 (x)x- 



(3.8) 



with a = 1/3. In figure ^3 we have plotted Too(h) / 'Nh a 
versus the scaling variable x = Nh e . The scaling hypoth- 
esis works fairly well. 



5 
z 

o 
H 




0.004 



0.002 




(Nh e ) 3 



Fig. 3.3. The scaling behaviour of the transition amp litud es 
Too (A', h) and T 20 (A, h) for small x [see Eqs. O) and (5t3)]. 



Moreover, the scaling function 

too(x)x — > foox, 



(3.9) 
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shows for small x-va lues a po wer behaviour, which is in 
agreement with ( |3.6|) and ( |3.7[ ). Finally, let us turn to the 
transition amplitude TW/i), which also vanishes in the 
limit h — > 0. In figure |3.3| we have plotted 



u T 20 (N,h) 
u; 2 o{N, h) 



(3.10) 



versus the scaling variable a; 3 . Here we have assumed that 
T20 and L02o/h scale in the sa me m anner. This is indeed the 
case as can be seen in figure 3.3. Moreover, in the small x 



r egio n the sca ling function behaves as it is expected from 
(I]!) and Q. 

In summary we observe that the excitation energies - 
with respect to the ground state - scale with h e (e ~ 2/3): 



w m o(N, h) = h e f2 m0 (x). 



(3.11) 



As before, we have neglected logarithmic corrections due 
to marginal operators. These corrections were deduced by 
means of conformal field theory (p3|. The transition am- 
plitudes scale with h a N (a ~ 1/3) 



T m0 (N, h) = Nh a mQ (x). 



(3.12) 



The small x-behaviour of the scaling functio ns f2 m p(x) 
and m o(x) is fixed by the initial conditions ( 2.13] ) and 
(pi: 



f2 m o{x) = a mQ x 1 [l + e m0 (x)], m = 0,1,2,... (3.13) 



According to (2.10) and (2.13) the small x-behaviour of 
different O m o{x) is different for m even and m odd, re- 
spectively: 



2 /m0(^) 



^° l J " I b m0 x-^[l + f m0 (x)} : (-1) 



(-1)" 1 = 1 



(3.14) 



4 The evolution equation in the scaling limit 



The critic al exponents a = 1/3 and e = 2/3 in Eqs. ( |3.1l| ) 
and (l3~ll ) as well as the small x -behavio ur 6 = 3 of the 
scaling functions shown in figures p||3~2j and pT3| will now 
be derived from t he e volution equations (2.2) and ( |2.3| ) in 
the scalin g lim it (3.3 ). We assume that a scaling ansatz of 
the type (|3~TT|) and |u|) 



uj mn (N,h) = /i e J? m „(x), 
T mn (N,h) = Nh°O mn {x) 



with the scaling variable 



x = Nh e , 



(4.1) 
(4.2) 



(4.3) 



holds for all the transition energies and amplitudes. Fur- 
thermore, let us assume that the small x-behaviour of the 
scaling functions: 



^mn(^) * ®mn% [1 



(4.4) 



„ / \ x^O J ClrnnX 2 fmn( x ) '■ (— l)" +m — 1 

mn(X) ~~ * 1 ^WT^l + / mn (.x)] : (-1)"+™ = -1 ' 

(4.5) 



follows from the initial conditions: 



and 



lim Nuj mn (N, 0) = a Tl 



T mn (N,Q) N -^F b mn N K 



(4.6) 



(4.7) 



where b mn — for (— l) m +" = +1. The exponents n,a 
and e are then related: 



k = 1 



(4.8) 



Again b m o is given by the initial condition (2.14). The 
large x-behaviour of the scaling function for the first ex- 
citation is of special interest. If 



lim f2 w (x) = fl\ 



(3.15) 



is finite, Eq. ( fc.ll ) tells us that there is a gap in the ther- 
modynamic limit, which opens with h e , e ~ 2/3. We have 
analysed our finite system data with the BST Q algo- 
rithm and found: 



10 



(3.16) 



A similar statement can be derived from Eq. ( |3.12 ) for the 
staggered magnetization, which can be expressed in terms 
of the transition matrix elements T m o(N, h): 



J(h) = J2\ T mo(N,h)f 



(3.17) 



Now we discuss the evolution equations fl2.4| ) and ( |2.5j ) 
in the sc alin g limit (|3.3|) , in conjunction with the scaling 
ansatz [(hi) and ( |4.2[) 1. It is convenient to inspect the 
derivative 



D 1 = N- 2 h-^+ 2 —u; mn (N,h). 



(4.9) 



Then the left- and right-hand sides of (2.4) acquire the 
following form: 



Di = e(e-l)x 2 n mn (x) 

+e(3e— l)x ; h e 



Dx = 2h 2 ^+ 2 



dx dx 2 

®LW , \- 2 n {x) 



(4.10) 



E 



E 



(4.11) 
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A similar calculation for the derivative: 

D 2 = h-^ In [T mn (N, h)u mn {N, h)] (4.12) 
yields for the left and right-hand sides of Eq. (pi|): 

D 2 = (<r + e) + ex-^ ln[0 mn {x)Q mn {x)], (4.13) 



Do 



x E 



®mn (^) 



l 



l 



fllm{x) Ol n (x) 



(4.14) 



The corres pond ing left-hand sides and right-hand sides in 
(4.9) and (4.12) scale in the same manner provided that 



2e- 1. 



(4.15) 



In this case (4T) and (12) lead to a system of differential 
equations for t he sca ling f unct ions f2 mn {x) and & mn (x). 
Combining ( 4.15| ) and (4.8) we see that the exponents 



6 = T~h? G = TT^' (416) 

are fixed by the exponent k in the initial conditions fl4.7| ) . 

The small a:- beh aviou r is governed by the initial con- 
ditions ( 2.14 ) and ( 2.13 ). They tel l us that the leading 
behaviour of the right-hand side of ( 4.11 ) arise from exci- 
tations with (-l)'+ m = -1 and (-1)'+" = -1 



x^o fii m {x) a im 



l-2er/e 



(4.17) 



Therefore, the terms in ( [4.11 ) should be proportional to 
Assuming that the small ^-behaviour of functions 
^mn(x) can be described by a power law ansatz: 



it follows from ( 4.10| ) that 

D\ — Q"mn&mnX^ 

is also proportional to x~ 3+2 / e if 

2 



(4.18) 
(4.19) 

(4.20) 



Therefore, the first coefficient e mn in the scaling function 
e mn (i) is completely fixed by the initial values ai n , bi m : 



2e mn = 2 (1 - (-1) 
1 



b 2 

m+n\ ^mn 
2 



(IT, 



E 



h 2 

ai m 



0-ln 



(4.21) 



Next we study the small x-behaviour of (4.12). We will 
discuss the two cases (~l) m + n = 1 and i^T) m + n = -1 
separately. 



4.1 = 1 



The right-hand side of (4.14) is governed by the excitations 
with (-l) l+m = -1 and [-!)'+" = -1: 



Omi{x)Oin(x) 

®mn (^) 
a;-»0 b m ibi rl 



1 



film{x) f2l„(x) 



1 1 



„2/e 



fmn (^) 



(4.22) 



For the left-hand side of (4.13) we get: 



3e — 1 + ex— \n[0 mn (x)Q mn (x)] 
ax 

^5 -1 + ex-^- bxf mn (x). 
ax 



Consistency of (4.22) and (4.23) is achieved if 



(4.23) 
(4.24) 



with <fif = <p e = 2/e. The coefficient f mn in (4.24) is com- 
pletely fixed by the initial conditions: 



f„ 



E 

l^m,n 
bmlhr, 



1 - (-iy+ m 1 - (-1) 



1 1 



l+n 



amn \ ai m ai n 



(4.25) 



4.2 (-1)™+™ = -1 

The contributions in ( 4.14| ) vanish: 
O m i(x)Oin(x) 



®mn (^) 



i 



Qlm{x) Oln{x) 



1 + {-l) m + l 


bi n a m i 


' 1 

f 


1 


2 


brim 


_ain 


aim 


l + (-l) n+l 


bmiain 


1 


1 


2 


bran 


_a m i 


a n i 



fml{x) 
fln(x), 

(4.26) 

since the scaling function fi m vanishes fo r the excitations 
with (— l)'+ m = l acco rding to equation ( 4.24 ). The lead- 
ing behaviour of ( 4.13 ): 



3e — 1 + ex— ]xi[G mn (x)Q mn (x)] 



j— o 



(4.27) 



is governed by the scaling functions e m n{x) and f m n (x ) 
for t he exc itations with (— l) m+n = -1. Combining (F2£) 
and ( 4.27 ) we obtain: 



2(6 mn ~\~ fmn*) 



l + (_l)i+m blnaml 
9 hnl , 



^ bmiain 
1^ o Jh 



i i 

a m i ai, 

1 ; 1 

a m i ain 



(4.28) 
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In summary we conclude: 
1. For consistency, all the scaling functions: 



.(*) 



. M) r 4>i 



(4.29) 
(4.30) 



have the same small x-behaviour with exponent <f> = 
2/e. 

Equations ( j4.21[) , (|4.25| ) and ( |4.2§| ) enable us to express 
the first coefficients e mn , f mn in terms of the initial 
values ai m , km- 

Having determined the lowest order in x of the scaling 
functions, we can proceed and compu te th e next orde r 
x^ 1 from the differential equations (4J3) and fl4.12| ). 
The exponent <j>i turns out to be <f>\ = 2<f) = 4/e for all 
the scaling functions. Moreover, the next order coeffi- 
cients emn, fmn can be expressed again in terms of the 
initial values a mn , b mn and the zeroth order coefficients 



^mn 7 J mn 



In 



The power ansatz ansatz (4.4) and (4.5) for the small 
a;-behaviour of the scal ing function is only consistent 
with a power ansatz (4.7) for the finite-size correc- 
tions of the initial values. Logarithmic corrections of 
the latter will induce logarithmic corrections in the 
small re-behaviour of the scaling functions. 



Finally let us discuss equations (4.21) and ( 4.28 ) for the 
lowest excitations eio and /io fro m th e gro und state. The 
numerical results shown in figures 3.1 and 3.2 tell us that 
the coefficients in the scaling functions for the higher exci- 
tations eo3 an d /03 are co nside rably smaller. If we neglect 
in equations (4.21) and (4.2J) all higher excitations we 
find 



eio 



,^10 

a 10 



and 

eio + ho — 0. 
A comparison of the slopes 

em = 0.0196(8), f w = -0.0177(8), 



(4.31) 



(4.32) 



(4.33) 



in figures |3.1| and |3.2| demonstrates that Eq. ( 4.32| ) is in- 
deed well satisfied. The numerical values of aio = 9.74 (8), 
610 = 0.95(3), es tim ated from the data of figure to- 
gether with Eq. (4.4) yield a value eio = 0.019. This is in 
agreement with Eq. ( 4.33] ) and might be seen as an indi- 
cation that the neglect of higher excitations is justified. 



5 The Heisenberg model in a 
dimerized field 

The general features developed for the ID spin-1/2 Heisen- 
berg model in a staggered field can be applied to all Hamil- 
tonians of the type 



H(<J) = H + SH.1 



(5.1) 



where H is the unperturbed Hamiltonian (1.1) and Hp 
denotes a perturbation operator. As a further interesting 
example we want to discuss here the dimer operator 



H, 



JV 

2 V^(— 1)™S„ ■ S n +i. 



(5.2) 



Hamiltonians of type ( p.l| ) with Hjj as (5.2) have been 
discussed in the framework of (organic) spin-Peierls mate- 
rials [p5|, where the interaction between linear antiferro- 
magnetic chains and the three-dimensional phonon system 
causes a dimerization of the lattice ]26],p7|,p8|. In case of 
the more recently discovered inorganic spin-Peierls mate- 
rial CuGe03, reference it turned out that in addition 
a next nearest neighbour coupling in the spin chain has 
to be considered pfl . We will not discuss the influence of 
such a term in this paper. 

The evolution equation for the energy eigenvalues E n (§) 
and the transition amplitudes: 



Tmn(N,S) = (!P m (a)|H D |!P n (a)> 



(5.3) 



are obtained from (2.4) and (|2.5| ) by a simple substitu- 
tion: h i-» S, E n (N,h) i-> E n (N,8), and T mn (N,h) i-> 
T mn (N,S). Note however, that the dimer operator ( j5^ ) 
conserves the total spin in contrast to the staggered spin 
operator. Therefore, the excitations from the ground state 
are singlet states as well. 

For the solution of the evolution equation we need the 
initial values for the excitation energies: 



lim Nujmn(N, S = 0) 



and for the transition amplitudes: 



,W0) 



for 



(— l) m +" = l 
for (-l) m+n 



-1. 



(5.4) 



(5.5) 
(5.6) 



T 
T 

-Lit 



Equation 

<5 = and the fact that the operator (5.2) changes the 



51) results from momentum conservation at 



momentum by 7t. Equation (5.6) is an ansatz for the finite- 
size dependence of the transition ampl itud es. Numerical 
data are shown for T w (N,0) in figure 5.1. They suggest 
an exponent k = 0.37(5), which deviates from the stan- 
dard values k = 1/2, predicted by conformal field theory 
[^9[|l^] . The discrepancy is due to strong logarithmic cor- 
rections, which are disregarded in the fit fl5.6| ). This leads 
to a shift of the exponent k. We expect that the standard 
value k = 1/2 will show up with increasing system size. 

Next w e st udy the finite-size scaling ansatz of the type 
(|rl) and (U|): 



uio(N, S) 
wio(AT,0) 
T W (N,S) 
T w (N,0) 



= 1 



eio (a) 
l + /io(x), 



(5.7) 
(5.8) 



for the lowest excitation energy oJio(iV, 8) and transition 
am plitu de Tio(iV, <5). Numerical results arc shown in fig- 
ure 5/2. Optimal scaling is achieved now with a scaling 
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variable 



N6 e 



0-72(1), 



(5.9) 



which is in good agreement with the value 

e= + ~ 0.73(5), (5.10) 
which follows from the exp onent k — 0. 37(5) in the initial 



value shown in figure 5.1. This value (5.10) is compati- 
ble with previous results extracted from a finite-size scal- 
ing analysis combined with renormalization group analysis 




0.04 0.08 
1/N 



0.12 




12 16 20 24 
N l/2 



Fig. 5.1. Finite-size dependence of the ground state en- 
ergy (a) and the corresponding transition amplitude (b) of 
Hamiltonian (fy}, for system sizes N = 8, 10, . . . ,24 (•). The 



solid lines represent fits to the asymptotic behaviour (2.13) 
and (Ell). 



The exponent <fi in the small x-behaviour of the s calin g 
functions eio(x) and fio(x) is in agreement with Eq. ( pMOj ): 



4> = 2/e = 2.74(5). 
The coefficients eio and fio in front: 

eio = 0.0346(8), / 10 = -0.0314(8), 



(5.11) 



(5.12) 



obey the approximate relati on (4.32). It is remarkable to 
note, that the relation ( |5.10 ) between the exponents e and 
k holds inspite of the fact that the true large N behaviour 
(k = 1/2, e = 2/3) is not yet visible in the finite system 
results. Affleck and Bonner [M have discussed the im- 
plications of logarithmic corrections to critical exponents. 
They have comprised those corrections into effective scal- 
ing dimensions i.e. critical exponents and found e = 0.78 
for chain length N = 20. 

In the same way, as we have a nalys ed the large x- 
behaviour of the scaling function ( |3.15 ) in a staggered 
field, we have extrapolated the gap in the presence of a 
dimcrized field: 



I/O 

2 



1.06 



1.04 



1.02 




0.96 



0.94 



f, (x)+l=-0.031(2)x 



(N5 E ) 3 



(N5 E ) 3 



Fig. 5.2. The ratios (5.7) and (5.8) for eio(x) and fio(x) versus 
{N5 e ) 3 with e = 0.73. The solid lines show the fits to the small 
a;-behaviour. 

6 Conclusions and perspectives 



In this paper we have developed a general procedure to 
investigate the behaviour of antiferromagnetic Heisenberg 
chains, which are weakly perturbed by an external field. 
As examples we considered here a staggered field and a 
dimerized field, respectively. 

The question of interest is: What happens with the 
excitation spectrum and transition amplitudes under the 
influence of the perturbation? The procedure used in this 
paper is based on the observation that excitation ener- 
gies and transition amplitudes for the perturb ation op- 
erator satisfy a system of differential equations (2.4) and 
(2.5), which describe the evolution in the stre ngth h (or 6 ) 
of th e ext ernal field. The initial conditions ( 2.1C ), ( ^.13| ) 
and (233) are completely fixed by the excitation spec- 
trum and transition amplitudes of the unperturbed system 
(h = 0, 8 = ). We hav e shown in section || that a scaling 
ansatz [(4.1) and ( |4.2| )1 leads to a consiste nt s olution of 
the evolution equations in the scaling limi t (|3.3|) . The ex- 
ponents a and e in the scaling ansatz [(4.1), (]4.2D] and the 



scaling variable Nh e are fixed b y the finite- size behaviour 
of the initial conditions ( pTC| ), fl2.13| ) and fl2.14| ). 

We have also tested numerically the predictions of the 
scaling ansatz and found good agreement for the expo- 
nents as well as for the scaling functions for small values 
of the scaling variable. Therefore, th e be haviour of the 
perturbed system in the scaling limit (3.3) is well under- 
stood. 

To answer the question whether a gap opens in the 
presence of an external field (h fixed) in the thermody- 
namic limit ( N — > oo), demands the knowledge of the seal- 



n? = 6.0(1). 



(5.13) 



ing functions (4.4) for large values of the scaling variable 
In principle, this asymptotic beha viour fo llows from a solu - 
t ion o f the differential equations ( 4.10 ), ( 4.11 ) and ( 4.12 ), 
(4.13) for the sc alin g functions using the initial condi- 
tions (4.4) and (4.5). In practice, it might be easier to 
study this limit numerically coming from finite system re- 
sults. Our numerical data are consistent with a constant 



behaviour ( 3.15 ) of the gap scaling function. 
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The opening of a gap in a uniform field has been re- 
cently observed J35| in a neutron scattering experiment on 
copper benzoat Cu(C6D 5 COO)2 ■ 3D 2 0. The existence of 
field dependent and field independent low energy modes 
at wave vectors q = tt(1 — 2M) and q = tt, respectively, 
leads to the expectation, that the compound is adequately 
described by the Hamiltonian (1.2) in a uniform external 
field BS^(p — 0). This system is known to show up zero 
energy excitations at wave vectors q = 7r(l — 2M) and 
q = 7T in the longitudinal and transverse structure factors, 
respectively @,^||- 

The exponential fit to the temperature dependence of 
the specific heat data revealed, however, that there is a gap 
which opens with the field strength B as B e , e ~ 2/3. This 
means of course, th at th e compound copper benzoat can 
not be described by (1.2) with p = 0. Oshikawa and Affleck 
ppf argued that the local g-tensor for the Cu-ions gener- 
ates an effective staggered field of strength h[= h(B) <C B] 
perpendicular to the uniform field B. Therefore, one is 
lead to investigate the Hamiltonian: 



H{h) = H + /iS 3 (0) +BSi(tt). 



(6.1) 



The authors of reference |36J studied the model (6.1) for 
_B = 0, which is just the case we investigated in this paper. 
They came to the same conclusion concerning the opening 
of the gap with h € , e ~ 2/3, as we found here. The ques- 
tion remains, how the exponent e changes if the uniform 
field B is switched on. Moreover, the gap might evolve in 
a different fashion for the field dependent and field inde- 
pendent soft modes q = tt(1 — 2M) and q = tt. We will 
address these questions in a further publication §37j. 
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